Exact dynamical exchange-correlation kernel of a weakly inhomogeneous electron gas 
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The dynamical exchange-correlation kernel f xc of a non-uniform electron gas is an essential in- 
put for the time-dependent density functional theory of electronic systems. The long-wavelength 
behavior of this kernel is known to be of the form f xc = a/q 2 where q is the wave vector and a is a 
frequency-dependent coefficient. We show that in the limit of weak non-uniformity the coefficient a 
has a simple and exact expression in terms of the ground-state density and the frequency-dependent 
kernel of a uniform electron gas at the average density. We present an approximate evaluation of 
this expression for Si and discuss its implications for the theory of excitonic effects. 

PACS numbers: 31.15.ee 



Since its introduction in works of Runge, Gross, and 
Kohn 0, Q , the time-dependent density- functional the- 
ory (TDDFT) has evolved into a powerful tool of in- 
vestigation of systems ranging from isolated atoms to 
bulk solids. In the important linear-response regime, 
the key quantity of TDDFT is the dynamical exchange- 
correlation (xc) kernel f xc defined as the functional 
derivative 



/xc[no(r)](r,r',w) 



SV xc [n] (r,w) 



5n(r',u>) 



n=n (r) 



of the dynamical exchange and correlation potential V xc 
with respect to the dynamical electron density n, taken at 
the ground-state value no of the latter. With this defini- 
tion, the density-response function \ can be represented 
in operator notation as 

X (r, r', u) = { [1 - xks{C + UcT'xks} (r, r', w) , (1) 

where xks is the Kohn-Sham (KS) density-response 
function of independent electrons, C = e 2 /|r — r'| is the 
Coulomb interaction, and e is the absolute value of the 
electron charge. While the density-response function of 
non-interacting electrons XKS can be straightforwardly 
calculated in many cases of interest (e.g., for homoge- 
neous electron gases in three and two dimensions it is 
given by the analytical Lindhard 's [3] and Stern's 0] for- 
mulas, respectively), the construction of f xc , whose role 
is to account for dynamical many-body correlations, is 
not straightforward. 

As an instructive specific case, let us consider the exci- 
tonic effect [5| in a semiconductor, which would manifest 
itself as an enhancement of the imaginary part of x f° r 
frequencies close to the fundamental absorption edge. We 
neglect for a moment local-field effects and write down 
the diagonal elements of the density response in momen- 
tum space as of Eq. (Q]) 



XKs(<i, q,^j 



l-Wq J q,w)[ 4 ^+/*c(q J q,w)] 



(2) 



where 47re 2 /g 2 is the Fourier transform of the Coulomb 
interaction. On the one hand, the excitonic enhance- 
ment of x is a many-body effect and, therefore, it needs 
a nonzero f xc to be accounted for within TDDFT. On 
the other hand, because of the divergent Coulomb part 
47re 2 /g 2 in Eq. ([2]), any f xc (q,q,uj) that remained finite 
at q = would give no contribution in the long-wave 
limit q — > 0. This simple observation shows that in order 
to include the exciton, f xc {<\, q, <jS) must be divergent in 
the long-wave limit at least as strongly as the Coulomb 
term. And indeed, when the q~ 2 divergence has been 
introduced empirically in papers dealing with the opti- 
cal absorption spectrum of semiconductors 0, H, it 
has yielded a good TDDFT description of the excitonic 
effect. 

Clearly it would be highly desirable to have a first- 
principle theory of the small-q behavior of the xc kernel, 
rather than relying on empirical parametrizations. In 
this Letter we take a step in this direction. We first show 
that the asymptotic relation 



lini / xc (q, q,u>) 



e 2 a(w) 



(3) 



(we introduce the e 2 so that a is dimensionless) is a rig- 
orous consequence of exact sum rules for the current den- 
sity response function. Then, in the limit of weak non- 
uniformity we obtain a simple and exact expression for 
a(uj) in terms of the ground-state density and the dy- 
namical xc kernel of a homogeneous electron gas at the 
average density [20(. 

We start by noting that the local density approxima- 
tion (LDA) to f xc is unable to produce the divergence. 
Indeed, within LDA 



/ K (r,r', U )=tKr),^(r-r'), 



(4) 



where f xc (n,ui) is the long- wave limit of the xc kernel of 
the homogeneous electron gas of density n. The latter is 
known @, 0, M, EH to be finite, no divergence arising, 
therefore, in the Fourier transform of Eq. ((4]). 
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To obtain an accurate non-local f xc , we resort to the 
recently proposed general method [12[ derived from the 
time-dependent current-density functional theory (TD- 
CDFT) [f3| • This method is based on the exact relation 
that holds between the scalar density-response function 
X (density response to a scalar potential) and the ten- 
sor current-density-response function \ (current-density 
response to a vector potential): 



X(q, q', = — j q • x(q, q', u) • q'. 



(5) 



Both response functions are expressed in terms of the cor- 
responding Kohn-Sham response functions and xc kernels 
in the following manner: 

x -1 (q,q',w) = XKs(q, q', - / ac (q,q',w) - -^-<W 



(6) 



and 



x (q,q',w) = XKs(q>q'> w )-/*c(q,q',w) — 2~£q<W - 

(7) 

where £q,ij = qiqj/q 2 , i and j are cartesian indices. 

Equations ([5][7|) establish a connection between f xc and 
its tensor counterpart f xc . The usefulness of this connec- 
tion stems from the fact that the tensor quantities xks 
and f xc satisfy a broader set of exact sum rules than the 
corresponding scalar quantities. These sum rules were 



derived in Ref. [14j. Specializing to the case of a peri- 
odic systems, the two most important sum rules for our 
purposes are 



XKS,ij(G,0,u>) 

1 V 



mc 

-^2 E XKS,i,k{G,G',u>)G , k G , j V K s{G') , (8) 

G',k 



■n (G)5, 



and 



V f xc ,ij(G, G',w)no(G0 = GjGjX^ c (G) , (9) 



G' 



where G are reciprocal lattice vectors. These sum rules 
connect three different types of components of, say, 
Xks{G, G', u>): the (0, 0) component, the (0, G ^ 0) and 
(G 7^ 0, 0) components, and the (G /0,G'/ 0) compo- 
nents. 

Let us further restrict our attention to the case of a 
weakly inhomogeneous system: |no(G)| <C no{0) = no, 
and \V K s{G)\ <C h 2 G 2 /2m for G ^ 0. Then it is eas- 
ily shown that the homogenous electron gas approxima- 
tion for the (G ^ 0, G' ^ 0) components completely de- 
termines the (0,G^0) and (G ^ 0,0) components to 
first order in no(G ^ 0), which in turn completely deter- 
mines the (0, 0) component to second order in no(G ^ 0). 
Thus, for xks we obtain 



^2 L G,ij Xks(G, w) + T Gtij Xks{G, w) 



XKS,ij(G ^0,0,u) = xksaj^-G,^ = ^[n {G)5 ij - L Gtij x h K L s {G,uj)V K s{G)] , 



0~GG' 
hL 



XKS,ij ■■{0,0, U}) 



en e 

°ij H 2~~2 

mc m z uj z c 



E G 2 L G ^\V KS (G)\ 2 [x h K L s {G,cj)-x h K L s {G,0)] 



(10) 



G^O 



to the zero-th, first, and second order in Vks{G), respectively. Here Xtcs ^^Xks are > respectively, the longitudinal 
and transverse KS density-response functions of the homogeneous electron gas of density no, and T Gj ij = $ij — L G ^j. 
Similarly, for f xc we have 



f xc , 3 {G^Q,G' 



G 2 [f^{G,cj)L G ^ + f^{G,u;)T G . l3 ] 6 GG , , 



fxdi(G 0, 0, w) = / XCiii (0, -G, to) ■■ 



cG 2 
euP-fiQ 



n (G) { {f X c (G, w) - f^{G, 0)}L G , ij + f^{G, lo)T g ^} 



/, c , tf (0,0, W ) = — ^ E G 2 \n Q (G)\ 2 {[/* C £ (G» ~ /^(G,0)] L G<ij + f^{G,cj)T G ^} , 



(11) 



G^O 



where /,£r and /^f are the longitudinal and transverse, 
respectively, xc kernels of the homogeneous electron gas 
of density no- 



The following steps, which involve repeated inversions 
of infinite matrices, rely on the mathematical fact that to 
find the (0, 0), (0, G ^ 0)/(G ^ 0, 0), and (G ^ 0, G V 
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(%KS> fxc) 
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f 

J xc 



FIG. 1: Scheme of the procedure for calculating the xc kernel 
/rcc starting from the expressions (|10f) and for Xks and 
f xc , respectively. 



0) elements of the inverse matrix to the second, first, and 
zeroth order in the inhomogeneity, respectively, it is suf- 
ficient to know the corresponding elements of the original 
matrix to the same accuracy, and_then the inversion can 
be performed in a closed form 



15] 



The complete procedure is schematically illustrated in 
Fig. Q. Starting from Eqs. (HUJ) and (HJ) for xks and 
fxc we (i) Invert Eqs. (fT0|) to get x~^ s ; (ii) Combine x~ks 
and f xc to get by virtue of Eq. 0; (iii) Invert 
to get x i (iv) U se Eq- © and its KS analogue to find the 
scalar response function x from x and Xks from Xks] ( v ) 
Invert x and xks to get and x~ks'i and > ( vi ) A PP!y 
Eq. |(7J) to find / xc . The final result of this procedure is 



nm/ xc (G ± 0,GV 0,w) - /i L (G, W )5 G& , 

km/ xc (G ^ 0,q,u/) = / xc (-q, -G,w) = [/^ L (G, W ) - ^ C L (G,0)] n„(G), 

q^o n q 

lim/ xc (q,q, W ) = ^ £ (G • q) 2 [/£f(G» - ^ C L (G,0)] |n (G)| 2 , 



(12) 
(13) 
(14) 



where q is the unit vector parallel to q. It should be noted 
at this point that the above expression for the scalar ker- 
nel /xc(q, q, u>) differs from what one would get by sim- 
ply taking the longitudinal component of /^^ (q, q, w), 
i.e. / xc (q, q, w) ^ ^ij '/</,-.,.. (q, q, u)qj ■ The im- 
plication is that the scalar xc potential (V xc ) of time- 
dependent DFT is not equivalent to the longitudinal com- 
ponent of the vector potential (A xc ) of time-dependent 
CDFT: rather, it should be constructed through the care- 
ful inversion procedure described above. A recent inter- 
esting attempt to construct V xc from A xc [l(| should be 
re-examined in the light of this result. 

From the result of the step (iv) for x and making use 
of the relation 



1 4ire 2 

— = 1 + lim — — x(q,q,w), 

e M (w) 9^0 q 2 



(15) 



we obtain a formula for the macroscopic dielectric func- 
tion of a crystal 



£m(w) = 1 



47re 2 n 



mw 1 



x £ |Vu(G)| 2 G 2 (q-G) 

G#0 



e hL {G,u) e hL (G,0) 



-(16) 



where Vq is the bare crystalline potential and 



JiL 



Aire 2 



Xks(<1' 



is the longitudinal dielectric function of the homogeneous 
electron liquid. Equation (|16p is in agreement with the 
Hopfield's formula for optical conductivity fl7j |. while in 
the RPA [f£ cL (G,u) = 0] it coincides with the corre- 
sponding result of Ref . |7j| 

Equations (|12j|14p are the main result of this paper. 
They replace the grossly inaccurate LDA formula 



lim f xc (G 



q, G' 



fLdG,u)6, 



GG'> 



(17) 



which does not contain the singularity in q. Identifying 
the (0, 0) element of the microscopic matrix of the xc 
kernel in Eqs. (fl"4|) as the averaged f xc , we see that f xc 
diverges for q — » as described by Eq. ([3]), wherein a(ui) 
is given by 



(G-q) 2 



[/^(G, W )-/^(G,0)]|n (G)| 2 . (18) 



G^0 



Notice that a(u>) = in the uniform limit and a(0) = 
up to second order in no(G ^ 0) [2lj ]. 

In order to calculate a(uj) from Eq. (fT5|) we need the 
Fourier amplitudes of the ground-state electron density 
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and the wave vector and frequency-dependent f££ of the 
homogeneous electron gas, evaluated at reciprocal lat- 
tice vectors. The first ingredient is straightforwardly ob- 
tained from standard electronic structure calculations. 
Unfortunately, the same cannot be said of the second in- 
gredient f^(q,uj), for which we do not have reliable ex- 
pressions. The best that can be done, at this time, is ei- 
ther to disregard the wave vector dependence, or to make 
use of the interpolation formula proposed in Ref. Uj£, 
which however fails to reproduce, at small q, what is 
presently believed to be the qualitatively correct form of 
the frequency dependence. In spite of these difficulties, it 
must be emphasized that the calculation of f^(q,Lu) is 
still a much simpler problem than the calculation of the 
dynamical xc kernel of the non-uniform system. Thus, 
our Eq. (|18[) does not simply express an unknown quan- 
tity in terms of another unknown quantity, but actually 
opens the way to systematic calculations of a based on 
the many-body theory of the homogeneous electron gas. 
Further, Eqs. (|12ll- (fT4|) for f xc offer a promising alterna- 
tive to the widespread practice of treating the dynamical 
exchange and correlations effects in the LDA. 

In Fig. [5] we plot a(u) from Eq. |T5J) vs frequency for 
crystalline silicon. The Fourier coefficients of the electron 
density were calculated with the code FHI98MD [l9( , and 
we approximated f^(q,uj) ~ f^(0,uj), taking the lat- 
ter from Ref. [HI In the range 0-22 eV, the real part of 
a(u>) is negative, changing sign for positive above 22 eV. 
It reaches its minimum of a ss -0.1 at lo ps 14 eV. In 
the range 3-5 eV of the main absorption in silicon, Re a 
changes from -0.01 to -0.03, which is an order of magni- 
tude smaller than the empirical value of a ~ -0.2 found as 
the best fit to the experimental spectrum in Ref. @. This 
large difference may simply indicate that the nearly free 
electron model, while being adequate for simple metals 
and even for semiconductors in the high-frequency regime 
[IH , is not sufficiently accurate for semiconductors at fre- 
quency lower than or comparable to the band gap (see 
also footnote [21]). Another probable source of discrep- 
ancy is that our approach is a pure TDDFT, whereas the 
value of a R* —0.2 was obtained in Refs.|6| and|7| with the 
use of self-energies incorporated in the Green's function 
via the GW approximation. 

In conclusion, within the nearly free electron ap- 
proximation, we have constructed the otherwise exact 
exchange-correlation kernel for time-dependent density- 
functional theory. This kernel is nonlocal in space, 
exhibiting the q~ 2 singularity in the reciprocal space. 
The strength of this singularity, which is frequency- 
dependent, has been directly related to the magnitude 
of the non-uniformity of the density of valence electrons, 
and this singularity disappears in the limiting case of the 
homogeneous electron liquid. We are proposing an im- 
provement over the conventional LDA scheme of includ- 
ing the dynamical exchange and correlation effects into 
ab initio calculations of the linear response of crystalline 
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FIG. 2: The frequency dependence of the real (upper panel) 
and image (lower panel) parts of the a coefficient in Eq. © 
for silicon calculated by Eq. (|18[) . 



solids which consistently accounts for the long-wave di- 
vergence in the exchange-correlation kernel. 
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